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PHOTOELASTIC  STRESS  AND  FRACTURE  ANALYSIS 
OF  TWO  NEUTRON  TUBE  DESIGNS 


INTRODUCTION 

A  neutron  tube  is  a  small  axisymmetric  component,  composed  of  inner  and  outer  metal  sleeves 
separated  by  a  glass  ceramic  insulator.  The  Sandia  Corporation  is  conducting  a  series  of  strength  tests 
on  various  neutron  tube  designs.  Test  results  indicate  that  premature  failures  are  occuring  at  or  near  the 
bond  between  the  inner  sleeve  and  the  insulator.  Examination  of  the  failed  specimens  indicated  that 
catastrophic  failure  is  proceeded  by  a  short  period  of  slow  crack  growth  at  the  interface.  In  order  to 
better  understand  the  mechanism  of  failure,  Sandia  has  conducted  an  extensive  series  of  analyses  on 
various  tube  designs. 

This  report  describes  stress  and  failure  analyses  on  two  candidate  tube  designs  by  the  Naval 
Research  Laboratory.  NRL  considered  one  loading,  a  uniform,  axial,  axisymmetric  load  acting  on  the 
protruding  end  of  the  inner  sleeve  and  reacted  at  the  opposite  end  of  the  outer  sleeve.  This  loading 
simulates  the  conditions  of  one  of  the  strength  tests  on  the  prototype  parts  conducted  by  Sandia.  The 
load  and  geometries  are  shown  in  Fig.  1. 

CONSIDERATION  OF  THE  METHOD  OF  ANALYSIS 

The  stress  field  in  the  neutron  tube  is  axisymmetric  since  both  geometry  and  load  are  axisym¬ 
metric.  The  materials  of  the  neutron  tube  are  the  ceramic  of  the  insulator  and  the  metal  of  the  inner 
and  outer  sleeves  which  has  a  Young's  modulus  approximately  four  times  that  of  the  ceramic.  The 
need  to  obtain  stresses  and  stress  intensity  factors  near  and  on  the  interface  of  the  two  dissimilar 
materials  is  the  most  difficult  feature  of  the  problem. 

Manuacript  submitted  Ausust  21, 1980. 
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SCHEMATIC  SHOWING  CROSS  SECTIONS  OF  BOTH  DESIGNS 


2D  MODEL  OF  M5N  2D  MODEL  OF  MI9N 

WITH  MACHINE  LOAD  WITH  AUTOCALIBRATION  SETUP 

Fig.  1  —  Geometry  end  loedinf  of  MSN  and  MI9N  neutron  tube  designs  and  models 
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A  number  of  methods  of  analysis  were  considered.  Classical  theoretical  methods  of  stress  analysis 
would  find  difficulty  in  representing  the  boundary  conditions.  Finite  element  methods  typically  give 
poor  results  on  and  near  the  interface.  There  is  no  practical  experimental  way  to  measure  the  actual 
stresses  on  the  interface  of  the  prototype.  Three-dimensional  bi-material  photoelastic  modeling  would 
be  costly  and  time  consuming.  Thus,  it  was  decided  to  conduct  an  approximate  analysis  using  two- 
dimensional  photoelastic  models  of  the  cross  section. 

The  rationale  for  approximating  the  axisymmetric  problem  with  a  plane-stress,  two-dimensional 
model  can  best  be  described  in  two  steps;  (1)  the  rationale  for  representing  the  cross  section  of  the 
axisymmetric  geometry  with  the  cross  section  of  the  corresponding  plane-strain  geometry,  and  (2) 
representing  the  plane  strain  cross  section  with  a  similar  plane  stress  model. 

First  Step 

The  solution  of  the  plane-strain  case  approaches  the  axisymmetric  solution  as  the  ratio  of  the 
inner  to  outer  radius  approaches  unity.  Clearly  in  the  neutron  tube  the  radius  ratio  is  not  near  unity. 
The  approximation  is  justified  for  the  neutron  tube  analysis  by  two  observations.  The  first  observation, 
based  on  the  principle  of  equilibrium,  is  that  most  of  the  stress  variation  between  the  plane  strain  case 
and  an  axisymmetric  case,  with  a  radius  ratio  far  less  than  unity,  is  due  to  the  gradual  increase  in  area 
of  a  representative  angular  (pie-shaped)  segment  of  the  axisymmetric  geometry  as  compared  to  a  uni¬ 
form  thickness  slice  of  the  plane  strain  geometry.  So  that,  if  the  same  force  is  applied  at  the  inner 
radius  of  the  axisymmetric  geometry  as  at  the  corresponding  position  on  the  plane  strain  geometry,  the 
in-plane  stresses  at  the  outer  radius  of  the  axisymmetric  geometry  will  be  less  than  at  the  corresponding 
area  of  the  plane  strain  geometry  by  approximately  the  ratio  of  the  inner  to  outer  radius.  A 
corresponding  reduction  will  occur  at  all  intermediate  points.  Thus  in-plane  stress  could  be  adjusted  to 
better  approximate  the  axisymmetric  case  by  a  simple  geometric  correction  factor.  The  second  observa¬ 
tion  is  that  the  variation  between  the  axisymmetric  and  plane  strain  stress  fields  cited  above,  and  any 
other  secondary  differences,  will  be  the  same  for  the  various  design  changes,  since  the  radius  ratio  is 
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fixed  in  all  cases,  and  the  radial  location  of  failure  is  the  same  for  all  designs.  Thus  for  purposes  of 
comparison  of  designs  the  stress  proportion  will  be  the  same. 

Second  Step 

The  plane-strain-plane-stress  correspondence  is  a  well  accepted  principle  in  elastic  analysis  and 
states  that,  given  certain  conditions,  the  stress  fields  in  plane-strain  and  plane-stress  are  exactly  the 
same.  One  of  the  conditions  is  that  the  problem  be  of  a  "first-boundary  value"  type.  The  metal-to- 
ceramic  interfaces  of  the  cross-section  violate  this  condition.  Thus  the  plane  stress  solution  will  not  be 
exactly  the  same  as  the  plane  strain  solution.  This  is  especially  true  in  the  vicinity  of  the  bonded  inter¬ 
face.  Here  again,  as  in  the  first  step,  some  engineering  judgment  must  be  used  to  decide  if  the  model¬ 
ing  gives  a  valid  approximation. 

The  stress  analyses  have  been  confined  to  the  free  boundaries  and  so  are  considered  good  approxi¬ 
mations  of  the  plane  strain  case. 

The  fracture  mechanics  analyses  on  the  other  hand  seek  to  determine  the  K-factor  on  the  bonded 
interface.  Fracture  mechanics  methods  have  long  dealt  with  plane-stress-plane-strain  correspondence, 
since  any  crack  because  of  it  very  small  tip  radius  must  be  considered  to  have  a  plane  strain  region 
about  the  tip,  whilst  away  from  the  tip  many  problems  have  a  plane  stress  character.  Between  the  crack 
tip  and  the  plane  stress  region  there  is  a  transition  region  from  plane  strain  to  plane  stress.  (This  tran¬ 
sition  region  may  also  contain  nonlinear  or  plastic  stresses  and  strains.)  To  avoid  the  transition  zone, 
(and  also  non-linearity  and  plasticity),  stress  measurements  are  typically  made  at  least  one-third  of  the 
model  thickness  from  the  crack  tip.  This  procedure  was  followed  in  all  the  analyses  reported  here. 
Thus  all  measurements  were  made  far  enough  from  the  bonded  interface  to  obtain  valid  analysis. 

The  rigidity  of  the  inner  metal  sleeve  was  simulated  with  a  pair  of  metal  clamps  on  either  side  of 
the  2-D  photoelastic  model.  A  similar  arrangement  could  have  been  used  to  simulate  the  outer  sleeve. 
This  was  not  done  because  it  was  felt  that  this  variation  would  be  influential  only  near  the  outer  sleeve, 
and  would  have  only  slight  effect  over  most  of  the  stress  field,  and  not  be  a  factor  on  the  stresses  at  the 
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inner  sleeve.  This  assumption  was  confirmed  from  2-D  finite  element  analysis  of  the  various  condi¬ 
tions. 

EXPERIMENTS 

Two-dimensional  photoelastic  models  were  made  of  the  cross  sections  of  two  neutron  tube 
designs,  M5N  and  M19N.  The  model  of  the  MSN  cross  section  was  16  times  actual  size,  and  the 
M19N  model  was  10  times  actual.  The  portion  of  each  model  representing  the  inner  sleeve  was  fixed 
between  a  pair  of  rigid  plates.  Axial  load  was  applied  to  the  portion  of  the  model  representing  the  outer 
sleeve. 

Light  and  dark  field  isochromatic  patterns  were  photographed  at  various  load  levels  on  both 
models.  After  the  isochromatics  for  the  stress  analysis  were  obtained,  the  same  models  were  used  for 
fracture  analysis.  The  models  were  removed  from  the  rigid  clamps  and  a  cut  was  made  along  the  line 
of  the  inner-sleeve-  insulator  interface  to  simulate  a  crack  where  fractures  had  been  observed.  This  cut 
was  made  with  a  thin  saw  in  some  models  and  an  angled  milling  cutter  in  others.  The  models  were 
replaced  in  the  rigid  clamp  and  reloaded.  Light  and  dark  field  isochromatics  of  the  loaded  model  with 
the  simulated  crack  were  photographed  as  before,  and  in  addition  close-up  photographs  of  the  crack  tip 
region  were  taken.  This  process  was  repeated  a  number  of  times  for  each  model,  lengthening  the  simu¬ 
lated  crack  at  each  step. 

ANALYSIS  AND  RESULTS 

Stresses  are  obtained  from  the  isochromatic  patterns  from  the  stress-optic  law 

*\  -  2  rm„  -  nf/t  (1) 

where  tr\  and  are  principal  stresses  at  a  point,  rm,t  is  the  maximum  shear  stress  at  the  point,  it  is  the 

fringe  order  at  the  point,  /is  the  stress-optic  coefficient  (a  constant)  and  ris  the  model  thickness. 
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On  free  boundaries  of  the  model  (in  particular  the  free  surfaces  of  the  insulator)  one  of  the  prin¬ 
cipal  stresses  is  zero  and  this  equation  reduces  to 

<7  —  nf/l  (2) 

In  the  dark  field  isochromatic  pattern,  the  fringe  orders  are  integers  (n  —  0,  1,  2,  3  ...)  and  are 
easily  identified.  The  light  field  fringes  are  the  half  order  values  (n  —  0.5,  1.5,  2.5  ...).  Intermediate 
fringe  values  of  n  can  be  obtained  with  simple  photoelastic  techniques  to  an  accuracy  of  about  0.01 
fringe.  Dark  field  photoelastic  patterns  from  the  MSN  and  M19N  models  are  shown  in  Figs.  2  and  3 
respectively.  Also  shown  in  Figs.  2  and  3  are  corresponding  numerical  solutions,  obtained  with 
TOTAL,  an  NRL  developed,  graphics-oriented,  finite  element  library  [1]. 

To  solve  for  the  stress  it  is  customary  to  determine  the  value  of  the  coefficient  /  using  a  simple 
calibration  model  with  a  known  stress  field,  of  the  same  material  and  thickness  as  the  model.  The  cali¬ 
bration  model  is  customarily  loaded  and  /  determined.  In  some  of  the  models  this  approach  was 
modified  by  putting  the  calibration  model,  a  tensile  strip,  in  the  load  train  with  the  neutron  tube  model 
so  each  was  subjected  to  the  same  load.  The  load  was  applied  as  shown  in  Figure  1. 

The  same  equation  applied  to  the  free  boundary  also  applied  to  the  tensile  strip,  since  the 
transverse  stress  in  the  strip  is  zero 

<rf  “  ncf/t  (3) 

where  the  subscript  c  is  used  to  indicate  calibration. 

The  load  in  both  the  tensile  strip  and  in  the  model  is 

P  -  (<rc)(Ac)  -  nef(wc  tc)/tc  (4) 

where  Ac  is  the  cross-sectional  area  of  the  tensile  strip  and  wc  is  the  width  of  the  tensile  strip. 

The  expression  for  load  reduces  to  r 

P  —  nc/we  (5) 
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Fig.  2  *  Dark  field  isochromalic  pattern  from  model  of  MSN  cross*section.  Finite 
element  solution  is  shown  for  comparison. 
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Fig  3  -  Dark  field  isochromatic  paticrn  from  model  of  Ml 9N  cross-section.  Finite  element  solution 

is  shown  from  comparison. 
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Combining  Eq.  (5)  with  Eq.  (2)  eliminates  / 

a  (6) 

t  n,  wc 

The  ratio.  Pit,  is  the  load  per  unit  thickness,  letting  P  -  P/t,  Eq.  (6),  for  the  boundary  stress  can  be 
written 


P 


n 

nc 


±_ 

wc 


(7) 


The  terms  of  the  equation  have  the  dimension  of  the  reciprocal  of  length.  Multiplying  the  equa¬ 
tion  by  a  representative  length  on  the  model  gives  a  non-dimensional  equation  and  a  non-dimensional 
(normalized)  expression  for  stress.  The  length  of  the  inner-sleeve-insulator  interface  is  designated  L 
and 

£4  -  JL  -L-  (8) 

P  "r  »>c 

This  approach,  referred  to  as  auto-calibration,  permits  determination  of  the  normalized  stress  without 
knowledge  or  measurement  of  model  load,  stress-optic  coefficient,  model  thickness  and  calibration 
thickness.  However,  these  values  were  determined  as  secondary  checks. 


The  normalized  stress  is  obtained  from  four  measurements;  n,  nc,  L  and  wc.  wc  is  the  same  for 
all  tests,  2.000  inches.  L  is  6.400  inches  on  the  M5N  model  and  4.000  inches  on  the  M19N  model.  nc 
varied  with  load  and  values  from  0.95  to  2.75  were  used.  The  n  value,  of  course,  varied  with  model, 
load  and  position.  Values  of  up  to  n  —  31.0  were  recorded. 

Once  the  a-L/P  ratio  has  been  established  for  any  point  on  a  model,  it  is  valid  at  that  point  for 
any  load,  and  any  size  model,  including  the  actual  neutron  tubes  it  represents.  The  L  for  both  neutron 
tube  designs  is  10.16  mm  (0.40  inches).  For  a  given  load  per  unit  length  P,  the  stress  at  a  given  point 
is  obtained  by  multiplying  the  <tL/P  ratio  by  P/L. 

The  normalized  stresses  on  the  free  boundaries  of  the  M5N  and  M19N  insulator  are  given  in  Figs. 
4  and  5  respectively. 
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Fig.  4  —  Normalized  tangential  stresses  on  Tree  boundaries  of  MSN  neutron  tube 
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The  stress  intensity  factors  (SIF)  for  various  crack  lengths  were  obtained  from  the  same  models 
used  to  obtain  stresses  on  the  free  boundaries,  and  some  of  the  analysis  was  similar.  Overall  views  and 

close-ups  of  the  crack  tip  for  the  MSN  and  M19N  models  at  various  steps  of  cracking  are  shown  in 

Figs.  6  and  7  respectively.  For  the  fracture  analysis  a  more  general  form  of  Eq.  (8)  gives  the  principal 
stress  difference 

<gjl  ~  <r2)L  _  _n_  L_  (9) 

P  r>c  wf ' 

The  principal  stress  difference  is  in  turn  related  to  the  opening  mode  SIF,  Kt,  the  shear  mode 
SIF,  Kn,  and  the  far  field  stress,  <rox,  by  a  form  of  the  Irwin-Sneddon  equations.  The  Irwin-Sneddon 
equations,  which  approximate  the  stress  field  in  the  local  neighborhood  of  a  crack  tip  (r/l  «  1)  are 
usually  written  in  terms  of  the  Cartesian  stresses 

1  L  0  I,  •  0  .  3#1  „  .  0 (,  .  0  30  H 

I  AT/  cos  y  jl  -  sin  ysin  —  -  K„  sin  —  2  +  cos  ycos  —  II  -  <r„, 

1  L  0  (,  .  .  0  .  30 1  .  „  .  0  0  30 1 

<Ty  -  1^/  cos  y  11  +  sin  ysm  —I  +  K„  smycos  ycos  — 

1  |„  .  0  0  30  ,  ^  0  |,  .  0  .  30 11  /im 

Tyy  -  yj—  I K,  sin  ycos  ycos  —  +  K„  cos  y  II  -  sin  ysin  —  II  (10) 

where  x  and  y  are  Cartesian  coordinates,  and  r  and  0  are  polar  coordinates,  with  the  origin  in  both  coor¬ 
dinate  systems  defined  at  the  crack  tip.  The  crack  has  length  /  and  is  located  on  the  negative  x-axis 
(0  -  180°). 

The  principal  stress  difference  is  related  to  the  Cartesian  stresses  by: 

(ff|  -  ej)J  ■  (<r,  -  e,)2  +  (2  r^)1.  (11) 


12 


NRL  memorandum  report  43^3 


****  lengths 


r 


PARKS  AND  SANFORD 


NRL  MEMORANDUM  REPORT  4363 


Substituting  Eq.  (10)  in  Eq.  (11)  and  normalizing  as  before  gives 


(<ri  —  <rt)I  :  | 

P  |  “  2wJ(r/I) 


K,y/vL  .  K„JirL  )  (  I  . 

- = —  sin  0  +  2 - = - cos  9  +  JCH  v  -  sin  • 

P  P  P 


(12) 


+  2  sin  (d/2)  \K,  >in  0{l  +  2coa9) 

P  «V2 rjL  [  P 


(1  +  2  cos2  0  +  cos  0)  J  + 

((f  1  -  )L 

For  a  given  point,  (r,  g),  with  a  known  value  of - = - determined  from  Eq.  (1),  this  is  an  equa¬ 


tion  in  three  unknowns:  Kt  •JvT/P,  Klt  VwT / P  and  <r„L/P.  Using  data  from  just  three  points  on 
the  model,  the  three  corresponding  equations  could  be  solved  for  the  unknowns.  Using  just  three 
points  could  lead  to  significant  error,  by  selection  of  insensitive  points  and  the  usual  measurement 
errors  [2].  To  avoid  these  difficulties  20  diverse  points  were  selected  from  each  light  and  dark  field 
photograph  in  the  vicinity  of  the  crack  tip. 


The  20  points  generated  20  equations.  An  overdeterministic  method  using  least  squares  [2,3]  was 
used  to  obtain  the  three  variables  for  eight  different  crack  lengths  on  each  of  the  two  designs  MSN  and 
M19N.  Several  models  of  the  MSN  were  tested  with  differing  procedures.  The  variables  were  com¬ 
puted  with  all  20  points  in  each  field,  but  were  also  computed  with  a  number  of  subsets  of  the  20 
points,  and  with  the  origin  shifted,  to  probe  for  any  significant  errors  in  the  analysis. 


All  the  SIF  values  obtained  were  checked  by  putting  each  set  of  Kh  Kp  and  a „  values  into  Eqs. 
(9-11)  and  calculating  the  fringe  order  at  each  measurement  point.  The  standard  deviation  of  the 

difference  between  the  measured  fringes  and  the  calculated  fringes  was  used  to  estimate  error  in  the 

% 

analysis.  For  the  M19N  data  only  3  out  of  16  sets  had  a  standard  deviation  greater  than  half  of  a  fringe 
order.  All  the  saw-cut  MSN  models  showed  a  standard  deviation  of  less  than  half  a  fringe.  On  the 
mill-cut  MSN  model,  S  out  of  8  had  a  standard  deviation  of  less  than  half  of  a  fringe  order. 
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This  same  approach  was  used  to  generate  graphic  displays  of  fringe  patterns  for  a  number  of  sets 
of  Ki,  K,/  and  <rM.  One  such  set  is  shown  in  Fig.  8  with  the  corresponding  photoelastic  pattern. 
Visual  comparisons  of  the  fringe  patterns  confirmed  in  all  cases  that  the  determined  SIF  values  did 
indeed  represent  the  analyzed  crack. 

The  opening  mode  SIF  and  the  shear  mode  SIF,  as  well  as  the  resultant  stress  intensity  factor, 
K-y/K?  +  AT//,  are  plotted  in  dimensionless  form  in  Figs.  9  and  10  for  the  MSN  model  and  in  Figs. 
11  and  12  for  the  M19N  model. 

DISCUSSION  OF  THE  RESULTS 

Stress  Analysis 

The  photoelastic  patterns  in  Figs.  2  and  3  and  the  normalized  stress  plots  in  Figs.  4  and  S  all  indi¬ 
cate  that  overall  the  stresses  are  greater  in  the  MSN  than  in  the  M19N  cross  section.  However  in  the 
critical  area,  at  the  end  of  the  inner-sleeve-  insulator  interface  where  failure  begins,  the  stress  in  the 
MSN  cross  section  is  reduced  appreciably.  This  is  largely  due  to  the  insulator  lip  in  the  critical  area, 
characterized  by  the  large  obtuse  angle  between  the  free  surface  of  the  inner  sleeve  and  the  free  surface 
of  the  insulator.  The  fairing  of  this  insulator  lip  attenuates  the  stress  in  the  critical  area. 

The  M19N  geometry  creates  just  the  reverse  effect.  The  overall  "arch"  shape  is  more  efficient  in 
transmitting  the  forces  between  the  inner  and  outer  sleeve  and  thus  overall  the  stresses  are  less  than 
those  of  the  MSN  cross  section.  However,  the  "arch"  creates  an  acute  angle  at  the  end  of  the  inner- 
sleeve-insulator  interface  that  multiplies  rather  than  attenuates  the  overall  stress,  and  so  creates  higher 
stresses  in  the  critical  area  than  found  in  the  MSN  cross  section. 

The  stress  analysis  indicates  the  M19N  cross  section  is  more  likely  to  fail  at  the  interface  than  the 
MSN  cross  section. 
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Fig.  8  —  Dark  field  isochromatic  pattern  and  pattern  regenerated  from  stresa 
intensity  factors  from  the  M19N  mode  with  a  length  of  UL  -  0.16.  Note  this 
pattern  is  viewed  from  the  opposite  side  to  that  shown  in  Fig.  7  (mirror  image). 
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PROTOTYPE  CRACK  LENGTH 

Fig.  9  -  Normalized  opening  mode  end  sheer  mode  stress  intensity  factors  of  MSN  neutron  tube  design 
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Fig.  10  —  Normalized  resultant  stress  intensity  factor  of  MSN  neutron  lube  design 
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Fig.  12  *  Normalized  resultant  stress  intensity  factor  of  M19N  neutron  tube  design 


21 


PARKS  AND  SANFORD 


Failure  Analysis 

The  failure  analysis  reinforces  the  stress  analysis.  The  MSN  stress  intensity  factors  in  the  critical 
area  are  lower  than  those  in  the  M19N.  Only  after  an  appreciable  crack  length  is  generated  in  the  M5N 
geometry,  so  that  the  crack  undercuts  the  insulation  lip,  do  the  M5N  SIF  values  approach  the  M19N 
values.  The  failure  analysis  also  indicates  that  the  M19N  cross  section  is  more  likely  to  fail  at  the  inter¬ 
face  than  the  MSN  cross  section. 

Factors  Beyond  the  Scope  of  This  Analysis 

Certain  factors  outside  the  scope  of  this  analysis  may  have  a  bearing  on  the  results  described  and 
so  are  included  here  for  completeness. 

(1)  Interface  failure  vs  material  failure 

The  stress  analysis  of  the  M5N  cross  section  showed  maximum  stresses  on  the  insulator  halfway 
between  the  inner  and  outer  sleeves.  At  first  glance  this  would  suggest  a  failure  away  from  the  inter¬ 
face.  However  failure  on  an  interface  between  two  materials  can  be  very  different  in  magnitude  from 
the  simple  failure  of  either  material.  The  following  comment  suggests  some  reasons  why. 

(2)  Interface  treatment 

The  type  of  bonding  of  the  inner  sleeve  and  the  insulator  has  a  great  influence  on  interface 
failure.  Factors  in  bonding  that  would  influence  future  include;  local  alteration  of  material  properties 
or  geometry,  addition  of  a  third  material,  or  introduction  of  residual  stresses.  If  these  factors  differ 
between  the  M5N  and  M19N  assemblies,  then  they  may  produce  results  that  run  counter  to  the  ana¬ 


lyses  presented  here. 
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